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Introduction 

The authors have tried to analyze procedures of the construction of differen- 
tial equations that are employed for modeling of macroeconomic processes. The 
results prove to be rather unexpected. Thus, the derivation of the differential 
equation of Harrod's model is based on a linear relation between capital and 
income. As a result, there arises a contradiction in terms of dimension that 
is rooted in incorrect treatment of the fundamental notion of the infinitesimal 
quantity. One can overcome this contradiction by relating capital to the integral 
of income over a corresponding time interval. However, in this case, the solution 
is by no means an exponential growth but a much more realistic relation that 
reflects, in particular, objective finiteness of the prognostic period. 

An analysis of the models of Harrod-Domar, Phillips, as well as of other 
models (see the well-known treatise by R. Allen), leads us to the conclusion 
that analogous deficiencies are, in principle, inherent in these models too. In 
general, the refraction in the sphere of economic dynamics of the methodology 
of the construction of mathematical models borrowed from the field of natural 
sciences, such as dynamics, electrodynamics, etc., proves to be absolutely un- 
justified. As a matter of fact, differential equations adequate to these models 
follow naturally from the consideration of an infinitesimal element. However, 
as regards the problems of economics, such an approach is objectively senseless. 
Nevertheless, economics, in its turn, has intrinsic advantages from the point of 
view of possibilities of mathematical modeling, which is embodied in the no- 
tion of balance. As we will show, there exist formal means to reduce Leontief 's 
model of "expenses-output" in its canonical interpretation to a system of linear 
differential equation (of, generally speaking, arbitrary order with respect to the 
derivatives). 

At the same time, the scantiness of the arsenal of the means of linear theory 
that are used in representative modeling of macroeconomic processes is almost 
universally recognized nowadays. In this regard, we will characterize briefly 
those areas of systems analysis that are devoted to the construction of non- 
linear models that are adequate to a given "input-output" mapping. In what 
follows, we nonetheless note that Leontief's model in the differential form can 
be elementary reduced to a Fredholm integral equation of the second kind (with 
respect to a vector function), whose theory and algorithms of numerical real- 
ization are as constructive as possible. In the case, when the kernel of such an 
equation depends on a parameter, which is quite naturally interpreted in terms 
of the object sphere, the spectrum of its possible solutions becomes extremely 
wide. We think that the development of the theory of Fredholm integral equa- 
tions of the second kind, whose kernels contain parameters, and its application 
to the modeling of the processes of economic dynamics is rather promising. 

Note on literature references in the English version: The reader should be 
advised that all the references to page and section numbers appearing in the text 
are given according to the Russian editions of corresponding literature sources. 

Contact address: voronin@hncu.edu.ua 



Chapter 1 

Analysis of Harrod's model 



1.1 Computational relations 

Let us turn to Harrod's model of the development of the economy represented 
in the book by L. V. Kantorovich and A. B. Gorstko [1] (pp.160, 161). In 
the authors' view, "despite the simplicity - the model takes into account only 
one limited factor, i.e., the capital - it can be used for a rude approximate 
investigation into the laws of the growth of the economy" . In the latter context, 
we focus our attention on the above-mentioned exclusiveness of the capital. The 
employed notation are explained as follows: Y (t) is the national income; K (i) is 
the capital (the assets); C (t) is the volume of consumption; S (t) is the volume 
of accumulation, and / (t) is the investment. All the quantities are dimensional. 
For the analysis of the model, it is reasonable to reproduce the text [1] . 

"It is assumed that the economy is functioning in such a way that the fol- 
lowing relations are fulfilled: 

Yit) = Cit) + Sit) (1.1) 

(the national income is distributed between accumulation and consumption); 

5(t) = /(t) (1.2) 

(the accumulation is equal to the investment); 

S{t)=^lY{t), (1.3) 

where fi — const (the accumulation constitutes a constant fraction of the na- 
tional income); 

X(t) = /(t) (1.4) 

(the growth rate of the capital is equal to the investment, where K = dK/dt); 

K{t) = uY{t), (1.5) 



where v = const (the capital-to-national income ratio is a constant quantity, 
which is an observed empirical fact). 

From the equations that describe the model, it follows immediately: 

Yit)^^I{t)^^kit) = ^Y{t) (1.6) 

p p M 

or 

Yit)/Yit) = (,/,.. (1.7) 

Consequently, the growth rate of the national income is equal to fi/v, and, 
therefore, if one assumes that at the initial moment of time (at i = 0) the 
national income is equal to Yq, the law of its change with time has the following 
form: 

y (i) = yoe''*/^ (1.8) 

which usually stands in quite good agreement with practice" . 

In obtaining relations (|1.6p and differential equation (|1.7|) , a key role is played 

by the dependence 

t 

K (t) = J I {r^) dr^ (1.9) 

-T 

that is obvious. Here, the capital K (t) can be measured only monetary equiv- 
alent whose unit will be chosen to be the dollar ($); t is the time variable mea- 
sured, for definiteness, in seconds (s). Given that dt has the same dimension as 
t, the function / (t) represents the intensity of the investment flow measured as 
$/s. A corollary of (|1.9p is just the derivative (|1.4p : there is simply no other 
derivative in the model (|l.ip - (|1.5p . 

On the basis of (|l.ip - (jl.3p . we can conclude that Y{t), C (t) and S (t) are 
also intensities of the flows of the income, consumption, and accumulation ($/s), 
respectively. However, there arises a contradiction, because Y (t) in (|1.5p is the 
national income per one- year period measured in $. Analogously, t in (jl.Sp 
denotes the number of years (their fractions). Therefore, the exponent fit/iy is 
dimensionless. Note that we had to carry out an investigation into the dimension 
of the components of the model. For the reason that will be explained in what 
follows, one does not specify clearly the dimension of employed quantities in 
macroeconomics of mathematical orientation. However, this issue underlies the 
most important point of the whole consideration! 



1.2 Refinement on the content of the model 

As Y (t) is the intensity of the income, relation (II. 5p had to look as follows 

K{t)^v f Y{ri)dr], te[U,t,+t,], (1.10) 

ti 



where i* is the time interval of one year; v is the number of years during which 
such an income counterbalances the capital of the j-th year. 

From a mathematical point of view, such a specification is of no great impor- 
tance: what is actually important is the fact that i* is finite. Behind this fact, 
one can easily see a transition to discrete analysis and the impossibility of deriv- 
ing any differential equation! This major objective of the present investigation 
will be developed below. 

At the same time, one can express the capital K (t) via the intensity of the 
income Y (t) in a way that is different from what is done in (|1.10p : namely, based 
on (|1.2p and (|1.3p . one can employ an analog of (|1.9p . Consider relation (jl.Sp 
strictly adhering to the definition [T] : 'V = const (the capital-to-national income 
ratio is a constant quantity, which is an observed empirical fact)" . In this regard, 
let us draw attention to the informativeness of the coefficient v (in contrast, /x is 
purely primitive). The essence is that the "observed fact" characterizes a basic 
coefhcient :^* for the time interval i* during which the income is measured. For 
i» = 1 year it is reasonable to set v^ — 10. Accordingly, for a certain time 
interval At, the coefficient is given by 

v = v,.Ull\t. (1.11) 

At that, the flow Y (t) ($/s) creates, in the period of time Ai, the volume 
of the income ($) equal to 

t+At 

Y{r])dr] (1.12) 

that, by the definition [1], is proportional to the capital K {t + At) with the 
coefficient ly: see (|l.lip . Accordingly, 

t+At t+2At 

K{t + M)^u I Y (??) d7?; K{t + 2At) ^u j Y {rj) drj. 
t t 

After differentiation, 

kit + At)^!^ [Y {t + At)-Y (t)] ; K{t + 2At) = v[Y{t + 2At) - Y (t)] , 
and, as a result of subtraction, 

k{t + 2At) -k{t + At)=iy[Y{t + 2At) - F (i + At)] . 
Assuming that the interval At is sufficiently small, we use a Taylor expansion: 

k{t + At)=k{t)+k{t)At + ...; Y{t + At)=Y{t)+Y{t)At + ..., 

where we have preserved the terms containing the first power of Ai and 2At. 
We get: 

k{t) = i^Y{t). (1.13) 



Thus, the considered model is reduced to the set of relations (|l.ip - (|1.4p and 
(jl.lSp . From these relations, it follows 

k {t) = tiY it) = vY {t) , 

which just yields fi = v, or, in the case fj, j^ v, 

where c is a constant. As regards the derivative Y [t) that, as in (|1.6p . leads to 
a differential equation, it has treacherously escaped. In general, the solution to 
the considered problem by means of (|1.12p can only be trivial. 

In other words, in order to equate the function Y (t) (with the coefficient) 
to its derivative in (|1.7p obtaining the exponential growth p.Sp . it was abso- 
lutely necessary to resort to the surrogate (jl.Sp instead of (ll.lOp . However, 
the "derivative" can be understood. Having used (|1.12p . we have essentially 
obtained an analog of the interrelation between the capital and the flow (|1.9p 
without complementing the model at the level of a quantitative content. These 
arguments are of heuristic nature. 

1.3 Improvement of the initial model 

Let us act in a different way, namely, by considering the time interval from to 
t. A priori, we assume that At is small. The income during this period is equal 

to 

t 

Y (r;) dr], 



and, by analogy with (jl.lip . 

v^vj^/t] (1-14) 

however, the assumption that v — const in relation (jl.Sp should be rejected. 
Indeed, it contains a contradiction resulting from the fact that u is related to a 
year [T] [whose duration i* is related to (jl.lOp ]. In order to remain within the 
framework of continuous analysis, we just use the rule of proportion. Thus, 

t t 

K (t) = ly I Y {r]) df] ^ ^^ I Y (tj) dr]. (1.15) 



By the way, for i ^ 0, the coefficient v —> oo, which is quite reasonable. 
Simultaneously, equation p.lSp takes the form K {0} — i^*t*F(0), and it can 
be easily compared to (|1.5p and (jl.lOp with regard to dimension. From p.lSp . 

it follows that 

t 

Kit)^-^jY{r,)dv+'^Y{t), (1.16) 





and, obviously, K (t) -^ oo for t ^ 0. By (II. 4p . / (i) ^ oo as well. This 
peculiarity will be eliminated in what follows by the cancellation of t. 
From relations (|1.2p - (|1.4p . we get 

kit) = ^IY{t), 

which, combined with (|1.16p . yields the equation 

l^J^)tY{t) = ^tY{t), t>0. 

Making a change of the variable 

i=t/U, (1.17) 



dt dt dt 

we get 



*■('■'- lis'' <'>»■ " = ^- "■'«' 

The solution to this equation has the form 

^(*>7-^- (1-19) 

(1 - at) 

On this basis, the following conclusions can be drawn: 

- if /i = (the absence of investment) or i/* — > oo in p.lSp . then Y (t) — Yq, 
which is not unreasonable; 

- the time interval of a reasonable forecast is reflected, because for i = a^^ = 
v^l \i the income function (|1.19p becomes senseless; 

- generally speaking, such a forecast is inherent in a reliable model, which 
should be contrasted with an infinite growth of the function (II. 8^ : 

- under the interpretation that /i — 0.5 (consumption and accumulation 
equally share the income) and v^ ~ 10 years, the period of a conditionally 
reliable forecast is also equal to 10 years, if it is set equal to 0.5cr^^. This seems 
to be realistic. 

It should be emphasized that Yq in (|1.19p is the intensity of the income for 
t = ($/s) rather than an income per year as in (|1.8p . Note that the coefficient v 
(as the number of years) is interpreted in sections 1.1 and 1.2, 1.3 as dimensional 
and dimensionless, respectively. This fact is a consequence of the observed in 
section 1.1 contradiction between the interpretation of Y [t) as the intensity of 
the income, according to (|l.ip - (|1.3p . and the volume of the income, according 
to (fTSl) . 

Thus, based on (|1.14p . we have obtained a rather satisfactory result. It is 
stipulated by qualitative different dependence of the capital K (t) on the flows 
7 (i) and Y (t): see (II. 9p and (jl.lSp , respectively. Namely, there is a variable 
coefficient t^^ in p.lSp . 



1.4 Discrete essence of the initial model 

Under the interpretation of Y (t), C (t), S (t), and / (i) as flows ($/s), the model 
is represented by relations (|l.ip - (|1.4p and p.lO|) . Using the definite integral in 
its simplest interpretation, we can represent (jl.lOp as follows: 

ti+t, 

K{t)^u f Y{r])dT] = i^UY{t^), te[U,U + U], i^O,f,.... (f.20) 

A change of the variables, by (|1.17p . leads to the following relations: 

Y(i)^C{t)+S{i); S{t)^l{i); S (t) ^ fiY (i) ; (1.21) 

k{i)^tj{i); (1.22) 

K (i) = vUY (U) ; t e [U,U + 1] , *; = 0, 1, . . . , (1.23) 

where Y (t), C (i), S (i), and / (i) are intensities of the flows measured in $/s. 
There is an obvious contradiction related to the fact that the functions in 
p.2ip - (|1.23p cannot, on the one hand, be discrete and, on the other hand, 
represent intensities of continuous flows. However, it is impossible to satisfy 
(|1.23p continuously, within the framework of the model (|l.ip - (jl.5p : see section 
1.2. By analogy with p.23p . one is only left with the option to represent relations 
pT^ and pr?^ in a discrete form: 

Y (£,) = C {U) + S (£,) ; S {U) = I (U) ; S (£,) = fiY (£,) ; (1.24) 

k{u)^i{u), (1.25) 

where 

Y (U) = t^Y (U) ; C (i,) = t^C (i,) ; S (u) = U.S (u) ; / (u) = tj (u) , 

(1.26) 
i = 0, 1, . . . ,n 

have the dimension of capital, which is quite unambiguously said in [1]: see 
section 1.1. 

At the same time, under inaccurate treatment, such an approach stands in 
disagreement with the rules of the calculus of infinitesimal that form the basis 
of differential models. Suppose that investments in years i, and U + 1 are given 
by li and li+i, respectively. Then, the capital is 



K,^K,^i+U; K,+i=K,^ 



ij+i- 



However, the definition K [tij = Li^i — Ii in discrete analysis is illegitimate. 

Should the techniques of continuous analysis be formally extended to a trans- 
formation of relations (|1.23p - (|1.25p . then, as it is said about (|1.6p (see section 
1.1), it "follows immediately" that 

dY 

-^ {U) /Y {U) = m/^, (1-27) 



and, accordingly, 

Y (U) = yoe'^*'/^ (1.28) 

which should be exactly the interpretation of the solution l\1.8\i . 

Thus, all the functions in (| 1.2311 - (1 1.25|1 arc discontinuous at t — ti. Other- 
wise, the system of these relations is degenerate, and the solution is trivial. As 
a matter of fact, discreteness stipulates a combination of factors: Y (t) ^ Yq\ 
the interval t* in (|1.20p is finite. 

The function K (i) is also discontinuous, because, by (|1.9p . (I1.17P and (|1.26[) . 

K (i) = I i{f])d7j = Ko + Y,^ (t^) ' < * < ^""+1 ' (1-29) 

-T/t, *=° 

where 



-T/t, 

As a consequence, this function is not differentiable in a usual sense at t — ti. 

However, by the existence of (|1.28p . the discontinuities of K (t^ at i = ii 
in (|1.29p are finite. Namely, they are analogous to discontinuities of the func- 
tion Y (i) with the coefficient ly. Summarizing, we arrive at the conclusion 
that the function (|1.29[) can be understood only as a generalized function, and, 
accordingly [2\ (pp. 57-60), 

n 

ii:(£)=^/(t;) (5 (£-£,), (1.30) 

i=0 

where 5 {t) is Dirac's delta-function, such that 

{CXD 
— OO 

By analogy with pT^ and p3(I|) . in prTT]) . 
1 



r(£) 



v 



Ko + Y^H^^ 



i=0 



<t <t 



n+l, 



dY (i) 



dt V 



n 



i=0 



and the functions C* (t) , S (t) , and / (i) in p.24p and (|1.25p have the same 
structure: in other words, all of them are generalized functions. At that, the 
ordinary differential equation with constant coefficients (|1.27[) has, indeed, a 
solution of the form (II. 8|) in the class of generalized functions [3J (pp. 60, 61). 



Can one, based on these arguments, draw a conclusion that 

Y (i) - Y (£) , (1.31) 

where Y (i) is determined by expression (|1.28p . and that the model (|l.ip - (|1.5p 
objectively reflects the dynamics of the macroeconomic growth? 

1.5 Inadequacy of the initial model 

In the explanation to (jl.Sp . it is said: "the capital-to- national income ratio is 
a constant quantity" yj. However, this fact predetermines the solution of the 
problem. Indeed, for t = 0, 1, . . ., we have, respectively: Kq, Yq = Kq/v; 

Iq = KniJ'lv; Kx = Xo (1 + a) ; ^1 = i^\lv; i\ = Kxiijv; a. = [ijv; 

K2 = Ko[l + {l + a) a] = Ko{l + a + a^) ; Y2 = K2/V] h = K2^^/iy; . . . ; 

71 

K„ ^KoJ2 «'; ^n = K,Jiy; /„ = Kn^i/iy. (1.32) 

i=0 

As a < 1 (this follows from the content of the model, which is, mathematically, 
insignificant), using the formula for a decreasing geometrical progression, we 
get: 

l-a"+i ~ l-a"+i ~ l-(n/i^r+^ , , 

y„-ifo^^ T=^o— ^y^ , , n = 0,l,.... (1.33) 

Such a result contradicts the solution (|1.28p that, for t = tn = n, has the 
form 

r„ = foe"" = foe'^"/^ (1.34) 

because the equality of expressions (|1.33p and (|1.34p would mean that 

„„ l-a"+i , l-a"+i 

e"" = — , an = In (1.35) 

1 — a 1 — a 

for any n — 1,2,... and < a < 1, whereas the cases n = and a = are 
trivial. 

Using an analog of (jl.35p for t = i„+i, we arrive at the equation 

1 - q;"+i 

a = ln- (1.36) 

1 — a" 

that is satisfied by a = 1. However, in this case, (|1.32p does not represent any 
progression, and, instead of (|1.33p and (|1.35p . we get: 

Y^^{n + l)Yo; e" = n + l. 
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This equation cannot be satisfied for any n = 1, 2, . . ., and it becomes appar- 
ent that the values of Yn determined by formulas (|1.33[) and (|1.34[) are cardinally 
different. Here, one cannot assume any approximation of the exponential depen- 
dence (|1.34p by means of the rational representation (|1.33p : see [3] (pp. 19-24). 
As a matter of fact, the same conclusion can be easily drawn by substituting 
concrete values of n and a into (|1.36p . 

However, a contradiction of the considered model is obvious: 

- on the one hand, the function Y (£) determined by expression (|1.28p satisfies 
equation pT?f|) : 

- on the other hand, the same function in the form (|1.34p is rather remote 
from (jl.33p that is, in essence, a programmed solution of the same problem in 
the difference formulation. 

Certainly, an explanation is related to the fact that both the functions Y (i) 
from section 1.4 and, accordingly, Y (t) in ()1.7p and ()1.8p can only be general- 
ized functions. In this situation, the equality (ll.31|) is impossible in principle, 
because, in reality, Y (t) is a number appearing as a result of an operation on 
the function Y [tj. In other words, there exists only a correspondence of the 
form 

Y{t)^ JY{i^)^{7^)dj^, (1.37) 



where (f it) is an infinitely differentiable functions on the interval [O, £„] except 
for its boundaries: it is called finite (alias trial or support) function. 
For example, 

(f\ - / ^^'^ i'^i/in) , i e [O, in] ; 
^ w - I 0, > i > i„, 

ip [tn/i] = 1. Moreover, the derivative in (|1.4p is understood as follows "T (p. 
34): 

k (-q) (p{r))dri = - K (77) ip (77) dr]. 



As regards expressions (|1.8p and (|1.28p . they can be associated only with the 
function Y (i) from p.37p . which by no means facilitates evaluation of Y (i) at 
the points t ^ ti. "But how can one define the integral of the product of a 
generalized function and a trial function, if one cannot work with values of the 
function at separate points? The answer is simple: in this case, one should 
define the integral axiomatically rather than constructively" [5] (pp. 10, 11). 
According to an alternative theory 6J, generalized functions are introduced as 
specially defined limits of series of continuous functions. From the point of view 
of reanimating the model p.ip - p.5p . they are also useless. 

I. M. Gelfand and G. E. Shilov point out: "In the solution of concrete 
problems of mathematical physics, the delta- function (as well as other singular 
functions) appear, as a rule, only at Intermediate stages. In the final answer. 
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singular functions are either altogether absent or figure under the sign of in- 
tegration in a product with some sufficiently good function. Thus, there is no 
direct need to answer the question what a singular function is in itself. It is 
sufficient to answer the question what the integral of a product of a singular 
function and a 'good' function is. ... In other words, we relate any singular 
function to a functional that puts this singular function and any 'sufficiently 
good" function into correspondence with a certain number" [3] (pp. 14, 15). 

However, why should we at all mention generalized functions, if their tech- 
niques are useless in our case? The whole point is that, only for these functions, 
the transformations leading to equation (|1.7p and, accordingly, its solution are 
legitimate. The solution of (|1.7p . understood in the classical sense, is continuous: 
see, e.g., [7] (pp.28, 29). At the same time, relations (|l.l|) - (|1.5p are fundamen- 
tally irresolvable in the class of continuous functions. One can draw a conclusion 
that they are inadequate to the process of the macroeconomic growth that is 
the subject of the model! What is more, the following arguments prove to be 
completely unsuitable: 

- macroeconomics operates the scales of decades; 

- against this background, an income per year is only a small step; 

- the function (jl.Sp smooths out such steps, which reflects the dynamics of 
growth on the whole. 

However, we have here not only just "steps": as a matter of fact, they 
induce a cardinal change of the type of the equation inherent in the problem. 
The inadequacy of the problem is stipulated, in the first place, by relation (|1.5p 
that provides it with the property of discreteness. In this case, the function 
Y {t) in ()l.l|l - p.5|) embodies fundamentally different qualities: the intensity of 
the flow and its scale on a finite time interval. The role of relation ()1.5|) is rather 
pragmatic: namely, as the derivative (|1.4p is objective, this relation relates it 
in a proportional way to the function, which leads to the differential equation 
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Chapter 2 

Other models 

2.1 Harrod-Domar's model 

In the work by R. Allen 8J (pp. 75-78), this model and the resulting dynamics 
of the growth are represented in a dimensionless form: 

Y{t)^C it) + I (t) ; C (t) = (1 - /.) F it) ; / (i) = z.^^, (2.1) 

where 0</i<l;i^>0. The author emphasizes that Y (t), C (t), and J(i) 
are " functions of continuously changing time" . The solution of the problem is 
obtained in the form 

YiFj^Yoe"'/". (2.2) 



For greater clearness, it is reasonable to transform (|2.ip and ()2.2|) into di- 
mensional notation: 

t = ht; Y (t) = YnY (t) ; C (t) ^ CoC (t) ; I (t) = hi (t) , (2.3) 

where Yq, Cq, and Iq are the intensities of corresponding flows ($/s). 
We get: 

Y (t) = hC (t) + k2l {t) ; hC (t) = (l-fi)Y (t) ; (2.4) 

k2l (t) ^ lytoYo, (2.5) 

where fci — Yq/Cq; k2 = Ya/Ia, which leads to the differential equation 

i^toY (t) = fiY (t) , (2.6) 

whose solution is 

Y{t) = Yoet''/'''°. (2.7) 

However, let us draw attention to the fact that the confusion of double-faced 
character of the function Y (t) from the previous section has disappeared, and 
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the situation has become completely transparent. Indeed, all the functions in 
(|2.4p are intensities of the flows; one cannot raise any objections against the 
derivation of equation (|2.6p and its solution (|2.7p . 

Nevertheless, expression (|2.7p constitutes a clear proof of the inadequacy of 
this model. The main point is the dependence of Y (t) on the dimension of the 
parameter tg from (|2.3p that, by definition, is chosen arbitrarily. In other words, 
a mathematical model, under the condition of its reliability, cannot depend on 
the scale of time. 

Accordingly, relations (|2.4p and (|2.5p do not contain objective meaning. Gen- 
erally speaking, the first two of them are quite elementary, whereas (|2.5I) . the 
so-called accelerator, essentially introduces a derivative. On can assume that, 
in this case, the problem was solved the other way, namely, in relationship to 
obtaining differential relations of the form (|1.6p . 

Note that the model (j2.ip does not employ the function K (t), the capital, 
whose special role is pointed out at the beginning of section 1.1 with reference 
to pp. Indeed, the dependence (|1.9p . as well as its corollary (|1.4p . can be called 
fundamental. The substitution of / (t) from p.4p into (|2.5p yields: 

k2K it) = lytoYo, 

from which we get 

Kiit) = {iy/k2)toYiit), (2.8) 

where 

Ki (t) = K{t)- Ko; Yi (t) ^Y{t)~ Yo 

are, respectively, an increase of the capital during the period of time from and 
t and the intensity of the income at the moment t. 

How can the capital in p.Sp be equal to the product (with a certain co- 
efficient) of the intensity of the income at the final moment of the term of 
accumulation and an arbitrary period of time? However, relation (j2.8p is, prac- 
tically, an analogue of (jl.Sp if one assumes to — t^, = 1 year. Accordingly, there 
appears discreteness (see section 1.4) that in the same way leads to a conclusion 
that the considered model is inadequate (see section 1.5)! 

It should be noted that R. Allen himself does not touch on the issue of 

dimension. It can be firmly grasped only on the basis of numerical examples. 

In this sense, the situation with the model of section 1.1 was more complicated. 

By definition 5] (p. 193), 

-,^ dK(t) , , 

/(t-) = -Ai. (2.9) 

Accordingly, for K (t) = K (t) /Kq, by analogy with ([^ . to reduce ((^ to the 
form (|1.4p . it is necessary that io — Kq/Iq. In this case, however, instead of the 
coefficient v, we get ^i^^ in (jl.Sp . If only the variable t is dimensionless, instead 
of (|1.4p . we get a relation of the type (|1.22p . This list can be continued. 

As is pointed out by T. Puu [9J (p. 76), exactly R. Harrod is the author 
of the idea of the formulation of the macroeconomic model in continuous time 
by means of a differential equation (1948). At the same time, an analysis of 



14 



the basic work by P. Harrod [TU], also published in 1948, does not confirm this 
fact. Computational relations in this work are given in a discrete form. On 
the basis of the same approach the model of Harrod, as well as the models 
of Domar, of Solow, and of Samuelson and Hicks, are considered by A. Lusse 
[TI] (pp. 163-166, 172-179, 180-182). L. Stoleru is of the same opinion about 
Harrod's model [H] (pp. 272-277). On the whole, one can assume that it is 
not R. Harrod to whom dubious merit of inventing the above-discussed models 
should be attributed. 



2.2 Phillips' model and other models 

K (t) is also absent from the accelerator-multiplier model of Phillips considered 
by R. Allen. We restrict ourselves to the equations of the general solution [S] 
(pp. 81, 82): 



dl{t) 
dt 



iit)-.'^^'^ 



dt 
dY{t) 



; Z(t) = (l-/i)y(f) + /(t); (2.10) 



= -X[Y{t)^Zit)] 



dt 

where, in addition to the functions employed above, Z (t) is the intensity of 
cumulative demand. The coefficients are specified as follows: 

K is the rate of reaction, i.e., the inverse of a constant lag of investment; 

I' is the factor of the accelerator power; 

^ is the multiplier; 

A is the rate of the influence of the production (income) on the demand. 
(Note that, in contrast to the model [T] where in five relations we had only one 
derivative, here, in only three relations we have three derivatives.) 

In the variables (|2.3p , the system of equations (|2.10p takes the form 

to/ (t) = -K [/ (t) - (i/to/fe) Y (t)] ; hZ it) = il-fi)Y (t) + k^I (i) ; (2.11) 

toY (t) = -X [Y (t) ~ hZ (t)] , 

where Z (t) = ZqZ (i); ks = Yq/Zq. 

Changing the variables according to (|1.17p . we reduce the problem to the 
solution of the second-order ordinary differential equation 

Y{t) + aY + bY{t)^0, (2.12) 

where a — ai/p; b — bi/ p^ , 

fli = K + /iA — KvX] hi = Ki>X, p = to/t^. 

It has the form 

y(£) =ciePi* + C2eP=*, (2.13) 
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where ci and C2 are arbitrary constants; 



Pi, 2 



1 



a^ ± 



ap - 



4&1 



6 = 0. 



are the roots of the quadratic equation p^ 

As in section 2.1 the transformations of the flows are quite correct. However, 
the situation is completely analogous to that considered in section 2.1. Because 
of the presence of the parameter p — to/t^ in (|2.13p . the solution depends on 
the choice of the scale of time, which makes the model inadequate. Otherwise, 
with the help of p, we could a priori set the periods of oscillations. 

Using pTlj) . we express / (i) in ([2ll|) via K (t): 



toK (t) = -n \k it) - {vto/k2) Y (i)J ; fcgZ (t) = {l-p.)Y (t) + k^K (t) ; 

toY{t)^-X[Yit)-k3Z{t)], 
which yields 

toK (t) = -kK (t) + (Kiyto/k2) Y (t) ; toY (t) = -pXY (t) + Afcai^ (t) 

As a result of the substitution 

fc2 



Yit) 



toK{t) + K{l-iyX)K{t) 



KVpX 

and some simple transformations, we get 

tl'K [t) + to{K + pX- nvX) k (t) + KvXk (t) 



0. 



(2.14) 



Here, absolutely inappropriate dependence of the solution on the scale of 
time is obvious. However, let us put aside this cardinal issue and assume that 
io = i*, as a law of Nature. It is not difficult to find that the coefficients of 
equations (|2.12p and (|2.14p completely coincide. At the same time, on the basis 
of (|1.4p . equation (|2.14p is a second-order differential equation with respect to 
the investment I{t). 

Accordingly, the functions Y (t) and / (f ) may differ from each other only in 
the constants ci and C2 in the representation of the solutions of the type (|2.13p . 
To find them, we need initial conditions 

Y (0) = Fo; Y (0) = %■ I (0) = /o; / (0) = h 

that are rather ephemeral. Combined with abstract character of the coefficients 
of equations (|2.10p . this fact makes the model practically useless. 

Besides, the considered model ignores, in fact, the function K (t), i.e., the 
capital, whose exclusiveness as a unique "restricted factor" is emphasized by 
the authors of [l]: see section 1.1. Indeed, dynamic process, by virtue of objec- 
tive circumstances, should be related to a comparatively stable factor, because 
otherwise, figuratively speaking, a reference point is lost. 
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Obviously, even under the condition of a " law of Nature" , the model (|2.10p 
has no rehabilitation potential. In other words, it cannot be converted even into 
a conditionally adequate one that would be analogous to the transformations of 
Harrod-Domar's model considered in section 1.3. 

In light of the above, let us turn to the same model of Phillips in the inter- 
pretation of L. Bergstom [T3] (pp. 40, 41): 



C{t) = il-^,)Yit); ~-^X 



=. ,-, dK (t) 



^^j[,Y{F)-K{F)], (2.15) 

where Y {t) is a, factual net income or output; C {t) is factual consumption; K (t) 
is the volume of the capital; n, v, 7, and A are positive constants (/x < 1). 

It is explained: "As the capital represents the 'stock', whereas output is an 
'outflow', the quantity v is inversely proportional to a chosen unit of time. Thus, 
for example, when time in measured in months, v is 12 times greater than in 
the case when time is measured in years." 

The problem is reduced to the solution of the equation 

-^ + (7 + /iA - i^jX) -^ + fijXK (t) = 

[a comparison with (I2.12|) yields 7 = k]. However, the model contains the fun- 
damental drawback that has been considered in detail above. In other words, 
equation (|2.15p is nothing but an analogue of (jl.Sp . Taking into account (|2.3p . 
for K {t) = KqK (i), we again arrive at the association of the capital ($) with 
the intensity of the flow Y (t) ($/s). Nevertheless, as follows from the above- 
mentioned quotation, the author of [13] fully admits a possibility of differenti- 
ating discrete flows. 

Having worked out certain techniques, we can a priori identify the compo- 
nents of macroeconomic models that contain contradictions. Thus, a "simple" 
version of Goodwin's model 8J (pp. 193, 194) involves the equation 

K{i)^vY (f) -f at : 

this situation is quite analogous to (|2.15p . 

From the point of view of the diagnosis, an "early" version of Kaletsky's 
model [8] (pp. 199-201), practically, does not differ much from the above. It is 
based on the relation 

/(t) ^ a{l - c)Y {t) - kK {T) + e, 

where a and k are dimensionless constants, or, taking into account a lag, 

1 ,^___ „^ ^___^, A 



Y 



^^-jji^)'^^^^+^^-^^^'\ + Y 
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Its "later" version [5| (pp. 205, 206), based on the equation 

X - .^ dY(i) dK (t) 

B{t) = a{l-c)Y{t) + 1^2 -^ - k^^ + e, 

where V2 and k are constants, mathematically, has practically the same defects 
as ((2l0| . 

Phillips' multipHer model [8] (p. 79) 

z(t) = (i-M)^(i); ^--A[r(i-)-z(i-)] 

is directly subject to the criticism of section 2.1. Note that the model is designed 
in such a way that it is impossible to access the function K (i), i.e., the capital. 

2.3 Differential and difference models 

Of interest are the arguments of T. Puu: "Samuelson's invention (1939) of 
the business cycle that combines the principles of the interaction between a 
multiplier and of an accelerator, is, without any doubt, a tremendous event. The 
fact that a combination of such simple reasons, i.e., the buyers' expenditure of a 
certain part of their income on consumption and the producers' preservation of 
a fixed relation between the capital and the production volume induced cyclic 
changes, was simple, surprising and, at the same time, convincing. This model 
has, so to say, scientific elegance. By the way, it should be noted that its 
economic prerequisite was the macroeconomic approach of Keynes. 

The initial model had been proposed as a process with discrete time, that 
is as a difference equation; later, it was rather skillfully developed in detail by 
Hicks (1950). Harrod (1948) arrived at an idea to formulate this process in 
continuous time, as a differential equation. The obtained system, instead of 
generating cycles, caused a growth; nevertheless, he clearly realized that the 
process of development was like balance on the edge of the blade submerged in 
surrounding instability. 

This established a tradition for several decades. Business cycles were formed 
on the basis of difference equations, whereas development was formed with the 
help of differential equations. Today, we could say that Samuelson and Hicks 
chose by chance a second-order process, whereas Harrod decided in favour of 
a first-order process. We could also admit that a choice between discrete and 
continuous modelling does not change anything. Dynamic processes of any 
order can be formulated in the form of difference equations and in the form of 
differential equations. 

A choice of the type of the model (discrete as compared with continuous) 
can be regarded as a matter of pure convenience. For analytical purposes, 
when we desire to apply theorems from vast literature on differential equations, 
a continuous approach is preferable. However, when we want to apply this 
model to an experimental time series that is inevitably discrete, we have to use 
discretization" [9^ (pp. 76, 77). 
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A. Bergstrom's opinion about the freedom of choice of the model is anal- 
ogous: "One of the most important methodological problems of constructing 
economic models is the question what equations should be used to describe 
such models: differential or finite-difference equations. Although many individ- 
ual decisions are made in regular intervals of time (say, once a week or once 
a month), variables observed by the econometrist represent a result of many 
particular decisions made by different individuals at different points of time. 
Moreover, the intervals of the observation of most economic variables are con- 
siderably larger than the intervals between decision-making represented by these 
variables. These facts lead to the idea that variables of a typical economic model 
should be regarded as continuous functions of time and that such a model should 
be described by differential equations. ... 

One more argument in favor of representing models in the form of differential 
equations is that, even in the absence of continuous observations of economic 
variables, the predictable continuous trajectories of changes of these variables 
may prove to be of considerable value. Let us assume, for instance, that, in the 
company's view, the volume of sales of its products is closely related to the na- 
tional income of the country. Then, to forecast the sales, it is very useful to have 
a predicted continuous trajectory of a change of the national income, although 
measurements of this variable are carried out only once a year. A continuous 
model allows one to obtain such a prediction using discrete measurements of 
economic variables over the past period of time" [T^j (pp. 8-11). 

The position of J. Casti is alternative: "In discrete time, the dynamics of the 
system can be described by means of difference relations. The most important 
property of such a description is that it gives us an idea of the behavior of 
the system in a certain local neighborhood of the current state. At that, it is 
implicitly assumed that local information can be somehow 'unified', which allows 
us to understand global (in time and space) behavior of the system. Such an 
approach proved to be sufficiently justified for an analysis of many physical and 
technical problems. However, the possibility of its application in the case of 
less studied problems, especially of systems of social-economic nature, is by no 
means obvious" [13] (pp. 17, 18). 

In addition, in light of the material of sections 1.1.-2.2, we have to refer to 
the remark of R. Allen of somewhat ambiguous meaning: "It should be noted 
that mathematical economics belongs to applied mathematics: it embodies a 
union of mathematics and economics. Any in the least bit interesting results in 
the field of mathematical economics can be provided only by an economist that 
uses mathematical techniques" [5] (p. 19). 

2.4 Abstracted model 

However, the above-quoted arguments [13], as well as partly [9j , implicitly imply 
an approximation of the behavior of the considered system that, in its turn, is 
a solution to an appropriate differential equation. In other words, the mathe- 
matical model is not constructed in an immediate relationship to an objective 
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content of the process. This fact motivated the choice of the title of the section. 
As an example, we turn our attention to a simplified version of the model of 
long waves [TS| (pp. 84, 85): 

X (t) = -p [x {t) - qy (i)] ; y (i) - -r [y {t) - sz (t)] ■ z{t)^x{t)-y (t) , 

where x (t) is the rate of an increase in labor productivity; y (t) is the rate of 
an increase in capital endowment; z (t) is the rate of an increase in the profit 
rate; p, q, r, and s are structural coefficients that can perform the following 
functions: 

- adaptation of the model to the behavior of the observed system; 

- estimation of related factors of the reliability of the model; 

- forecasting of the behavior of the system (under the condition of adequacy 
of the model). 

In [15], it is pointed out: "The system of differential equations with positive 
and negative feedbacks employed by the authors describes a process analogous 
to a simple mechanical system such as, e.g., a pendulum, a spring, elastic con- 
structions with a damper, etc. However, there exists a considerable difference 
between economic and mechanical systems. In the simplest mechanics, elas- 
ticity and damping forces are usually independent of the system, whereas in 
economics all the factors are mutually interrelated" . 

In particular, the following is established. Regular cycles appear for s = —2 
and p = r. In the interval p e [0.10,0.12], periodic undamped oscillations 
appear whose duration is 50-60 years. Cycles of twenty years appear when p 
and r increase up to 0.34. A more complete version of the model was tested 
using the statistics of time series of the USA during the period from 1989 to 
1982. The following values of the "adaptation coefficients" have been obtained: 
p = 0.048 and r = 0.25, which corresponds to damped oscillations with the 
period of 53.7 years. 

The considerable difference between the coefficients p and r is objectively 
interpreted. However, the sensitivity of this model to a small change in the 
coefficients is apparent. From this it follows that the model, represented by a 
system of two first-order differential equations with constant coefficients, has 
rather limited the potential of the description of macroeconomic processes. 

2.5 Leontief s model 

In R. Allen's view, "the economist can learn much from the engineer: both the 
ways of using mathematical methods and the ability to pose technical problems" 
(p. 19). He can be well understood: indeed, equation (|2.12l) of Phillips' model 
contains four absolutely abstract coefficients k, v, /i, and A. In addition, we have 
"independent investment and consumption expenditures". They represent the 
free term of the equation and, accordingly, should be given. It goes without 
saying, that the fundamental defect caused by the dependence of the solution 
on a choice of the scale of time is also present: see sections 2.1 and 2.2. 
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However, equation (|2.12p is widely used in technics. Thus, it describes free 
oscillations of the mass suspended to a spring under the condition of viscous 
drag. All the parameters of such a system are concrete, they can be measured, 
and the same concerns the external force in the case of forced oscillation. The 
differential equation is rigorously derived on the basis of the laws of mechanics: 
see, e.g., [H] (pp. 43-49). 

So, what can be learnt from the engineer, if we aiming at developing an 
analogous approach to mathematical modelling in macroeconomics? It seems 
that anything like that is impossible by definition. In economics, there are no 
laws for idealized objects that can be put into correspondence with a material 
point. However, economics, in its turn, has great advantage over mechanics that 
is embodied in the equation of the balance of financial flows! 

Namely, Leontief 's model is brilliant both with respect to its simplicity and 
efficiency related to the development of computer technologies. For functions 
that depend discretely on time, this model was studied by M. Morishima. It was 
pointed out: "So-called Leontief 's model, initially static by its nature, is usually 
transformed to a dynamic one by the insertion of consumption and output lags, 
of the lifetime of the means of production, of accelerators, of the growth of final 
demand, of structural changes caused by technological innovations, etc." [T7] 
(P- 77). 

One can draw a conclusion that the possibilities of mathematical techniques 
of mechanics and economics, figuratively speaking, have different vectors. The 
main thing is their reasonable application rather than following the principles of 
blind imitation of approaches from those fields of knowledge that seem to be the 
most mathematized ones from a superficial point of view. By the way, although 
the engineer is aware of exact solutions for a bar, a plate, and other elements 
of this kind, he experiences considerable difficulties when designing unique con- 
structions. He has no restrictive criterion for the whole set of elements. The 
balance is absent! 

In light of the above, proceeding to the constructive part of our consideration, 
we consider, following T8] (pp. 40-45), the dynamic of a system defined as 

Xit) = F{Xit-T)), (2.16) 

where X (t) is a vector function of a set of flows Xi {t); Tm = max{ri} is a 
maximum delay (lag) during which a change of the state can be taken to be 
linear. 

The dynamics of the system is characterized by the derivative of the function 

■ . . dF . ^ dF , , , ^ 

^(t) = ^X(t-r) + -X(i-r). (2.17) 

From here, taking into account that r„j is small, we can get: 

^ W = — ^ [^ W -X{t-r)] + ^X {t-r), (2.18) 

Tm dX at 

or 

X{t) = Git,X{t),X{t-T)). (2.19) 
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At the same time, because of the presence in (|2.18p of the small factor r^ 
by the derivative, equation (|2.19p belongs to the class of singularly perturbed 
equations. Both the study of such equations and their numerical realization 
require the use of techniques from a rather special arsenal: see [12] (sections 2, 
7) and also [ID], [H] (section 1). 

In the next section, we present a method of the reduction of Leontief's 
balance model to a differential form which is alternative to the use of (|2.16p and 
dUTT]) that lead to (^1^ . However, if the function F (X) in ([^1^ is hnear, 
equation (|2.18p becomes meaningless by virtue of its triviality, and, accordingly, 
the derivation of a differential equation proves to be fundamentally impossible. 
In this regard, we note that the issue of limitations of a linear theory has been 
indirectly touched on also in section 2.4. This fact is related to the content of 
section 3.2. 
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Chapter 3 

Constructive arguments 

3.1 Leontief s differential model 

Consider a system of balance equations: 

n 

X, (t) = J2 «u2;j (i) + c, (t) , i = 1, 2, . . . , n; < i < io (3.1) 

that can be interpreted, e.g., as follows: Xi (t) are the flows of the volumes of 
output ($/s); aij is a part of the commodity i used in the production of the 
commodity j; Ci [t) > is the flow of an external demand for the commodity i; 
t is the variable of time (s). In what follows, the statement of the problem will 
be complemented. In matrix notation, the system of equations p.l|) takes the 
form 

X (i) = AX{t) + C{t), < t < to. (3.2) 

In this case, it is reasonable to assume that the matrix A satisfies the con- 
ditions of Metzler's theorem [17] (p. 34), and, accordingly, 

n 

^a,y < 1, i = l,2,...,n, 

where we have a strict inequality for at least one of the sums. Then, equation 
(|3.2p can be solved by the method of simple iterations [H] (pp. 120-122): 

X,+i (t) = AX, {t) + C{t), s = 0, 1, . . . , (3.3) 

where, in theory, the initial element can be chosen arbitrarily. 

However, we have the initial condition X (0) — Xq, and, accordingly, the 
iterations (|3.3p have, in reality, the form 

X (ts+i) = AX its) + C (t,) , (3.4) 
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where tg+i = tg + At, and there emerges the issue of choosing the step At that 
corresponds to the factors of changes in the vector function X (t). 

To approximate this function, one employs a Taylor series expansion and 
usually retains the first-order derivative [23j (pp. 32, 33). However, if we only 
want to capitalize on the idea of the outlined approach and, instead of 
consider the relation 



Xit + to) = AX{t) + Cit), 0<t< to, 



(3.5) 



where the time interval ig is sufficiently large, we will have to retain in the 
Taylor series 



x{t + to) = x {t) + toX (t) + Ulx (t) + hl'ji (t) 



(3.6) 



more terms. 

If, as an illustration, we restrict ourselves to three terms and use the notation 
A = E - B, where 



B 



I 1 - an ai2 

021 1 - 022 



\ a„i 



0.112 



1 



0.2n 



(3.7) 



and 



(X 
1 







Vo ... 1 / 
is a unit matrix, after substitution into (|3.5p . we obtain 

t^X (t) + 0.5^0^ {t) = -BX (t) + C (t) 

that under a change of the variable t — t/to becomes 

0.5X it) + X (t) + BX{F) = C{F), < f < 1. 



(3.8) 



(3.9) 



As a matter of fact, a derivation of this equation has been the aim of the pre- 
vious transformations. In the solution of this equation, one uses initial condition 
of the form 

X(0) = Xo, X(0)=Xo, (3.10) 

and, despite its matrix form, the procedure of numerical realization is quite 
standard [2^ (section 10). This reference also gives the solution by quadrature 
of the first-order differential equation with the elements of the matrix ()3.7p 
depending on time, i.e., for Uij = aij (t). In principle, this result can be easily 
extended to the case of the second order as well as of higher order: see, e.g., [7] 
(section 5). 
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However, one may ask whether we are resorting to a double standard here 
because the procedure of the reduction of (|3.8p and (|2.1ip to the forms p.9p and 
(|2.12p . respectively, is the same. Moreover, why the dependence of the solution 
on the time scale has not been pointed out in this case? This is a very important 
issue resulting from the use of the Taylor series (|3.6p . The number of retained 
terms in this series directly depends on the value of to. 

The difference lies in the fact how this dependence is established: by a trial- 
and-error method or by a profound analysis. As regards the consideration of 
sections 2.1 and 2.2, there was no criterion of this kind there. 

In the formulation (|3.ip . as the function C(i) is not given, the demand 
is prevailing and, in practice, the output may fail to meet it. To avoid this 
situation, we can multiply C (i) by a factor < a < 1 defined, e.g., by the 
condition 



f X{r])drj = X, 



where X^ characterizes the demand for internal output. To carry out analyt- 
ical studies for forecast purposes, one can complement the vector C {t) by a 
component that nonlinearly depends on X (i). 

3.2 Nonlinear model and an alternative 

"The author has arrived at a final conclusion that linear dynamic modelling 
yields very little. This is caused by the scarcity of the set of alternatives, i.e., 
either damped or explosive motion, that are related to linear models. Therefore, 
in the present study, we focus mainly on nonlinearity" [9 (p. 7). In what follows, 
we show that nonlinearity has its own alternatives. 

In our view, the work of A. van der Schaft 25j is of great interest to math- 
ematical modelling. It contains information about the orientation of system 
theory concerned with the realization of nonlinear dynamic models. Thus real- 
ization implies the construction of a system of equations 

i = / (x, u) , X (0) = xo, y = h (x, u) , (3.11) 

where x{t), y{t), and u{t) are vector functions, that is adequate to a given 
input-output mapping 

y{t)^F{u{T);0<T<t). (3.12) 

In other words, corresponding measurements can be carried out with a shift in 
time, which is rather important. All the quantities are dimensionless. 

As is pointed out by J. Casti, the problem of realization consists in the 
construction, as far as possible, of a compact model that stands in agreement 
with observed data [14]. 

Note that, for analytical studies of the nonlinear model (|3.1ip .. there exist 
efficient mathematical techniques dating back to the works of A. M. Lyapunov 
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and A. Poincare. The monograph [5S] is devoted to the adaptation of these 
techniques to the problems of macroeconomics. 

The author of [5S] has formulated an original part of his work in such a 
way: "We consider systems of smooth nonlinear differential and algebraic equa- 
tions where certain variables are singled out as 'external' ones. The problem 
of realization amounts to the substitution of an implicit high-order differential 
equation by explicit first-order differential equations and algebraic equations by 
means of a mapping in terms of external variables" . 

Thus, instead of (|3.12p . a relationship between the vectors at the input and 
output of the system, i.e., u{t) and y{t), respectively, is implied in the form 
of a nonlinear high-order differential equation. Certainly, with regard to the 
problem we are considering, such formulation is absolutely irrelevant. 

Indeed, in the course of consideration, we have been trying to derive a dif- 
ferential equation of at least the first order, and we have been discussing the 
correctness of the application of this procedure in the works on macroeconomics 
that are considered to be classic. At the same time, the mapping (|3.12p is rather 
organically related to the factors of the balance of financial flows, to the realia 
of corresponding measurements as well as to other factors of this kind. It may 
seem that one can be satisfied with this fact. Why have we then referred to the 
above quotation? 

The reason is that the text [25j implies an alternative way of the construction 
of macroeconomic models. Indeed, as a prerequisite of further transformations, 
here appears an input-output mapping that is represented by a system of linear 
differential equations 

^(1)^W=^(I)^(^)' ^'-''^ 

where D (s) and N (s) are polynomial matrices of "appropriate dimensions", 
subject to a number of loose requirements. This case is matched by the model 
(j3.1ip under a linear approximation: 

X ^ Ax + Bu; Y = Cx + Du, 

where A, B, and C are corresponding matrices. 

Consider, however, an ordinary differential equation, with constant coeffi- 
cients as in (|3.13p . of order n: 

z(") + a„_iz("-i) + ... + aiz^O, (3.14) 



where z(") = {d/dtf z{t). 
Using a formal notation 



we get 



t 



z^-Ht)^^{t). 



'^'^ = j^J^'-'^y'^'^^'^^^'^ + T.^^^'^j^: (3-15) 

*-i 
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Here, Ci are the constants of integration defined by conditions at the input- 
output of the system, that is, at i = and, for definiteness, at < = 1. 

The substitution of the function (|3.15p and of its derivatives into (|3.14p leads 
to a Fredhohn integral equation of the second kind in the canonical form: 



ip{t) = X k [t, ri) if (77) d77 + g (i) , t e [0, 1] , (3.16) 



where A is a parameter discussed below. As a matter of fact, we have described a 
classic method of the reduction of an ordinary differential equation to an integral 
equation that is given in numerous publications: see, e.g., [27] (pp. 31-33). 
Note that the use in (13.81) of the notation 



yields 

t 

X{t)^ I [t~vi)U [tj) dTJ + Cot + Cl, 



which, under conditions (|3.10[) . results in Leontief's model in the form of a 
Volterra integral equation of the second kind. 

However, in light of the present consideration, the following point is of pri- 
mary importance. The whole information about the system, i.e., about equation 
(|3.14p and conditions at i = and t = 1, is reflected by the kernel k {t, if) and 
the free term q (i) of equation (|3.16p . At the same time, as a result of the re- 
peated integration according to p.lSp . the kernel in (|3.16p represents, in reality, 
a rather particular case, i.e., a convolution 

k (t, ?/) = fc (t - 77) , 

and this convolution is of a very specific type. 

In this regard, one may ask a natural question: Based on a priori information, 
measurements, experiments as well as other factors including heuristics, why not 
construct the kernel k (t, 77) in such a way that the function tp (t) could possess 
the property to represent potentially realizable variants of the behavior of the 
considered system? 

Simultaneously, the kernel k (i, ?/) should contain parameters intended for the 
purposes of adaptation and correction of the model, as well as of forecasting. 
They can be associated with a part of conditions at t — 1 for the function z (t) 
and its derivatives from p.l4p that are subject to identification. Without any 
doubt, such a model contains elements of self-training. 

However, could a model based on the integral equation (|3.16p prove to be 
more efficient than (13. lip , that is, of the model posed in the form of the Cauchy 
problem that used to be a traditional orientation of the macroeconomic science? 

From this point of view, E. Goursat's constatation is of interest: "The so- 
lution of Volterra's integral equation is a broad generalization of the Cauchy 
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problem for a linear differential equation" f55^ (p. 16). At the same time, 
however, equations of Volterra's type represent a particular case of Frcdhom 
equations whose spectrum of solutions is more diverse. 

To answer the above-posed question, we consider the representativeness, in 
a mathematical and, accordingly, objective sense, of the solution of equation 
p.l6p that depends on the parameters of the model. In this regard, solutions of 
the exponential type that were present in sections 1 and 2 can be called trivial. 

3.3 Representativeness of the integral model 

Note that, for every set of the coefficients a^ from (|3.14p . the parameter A in 
equation p.l6p is a definite number. At the same time, to study the properties 
of the function tp [t) that satisfies (|3.16p , it is reasonable to consider A indefinite 
and to specify it when necessary. 

The solution to equation (|3.16p in the case of a sufficiently small absolute 
value of the parameter A can be expressed via the resolvent: 



ip{t)=q{t) + \JH{t,r^,\)q{ri) 



df], 



under the condition A ^ Ai , where A^ are the characteristic numbers (alias points 
of the spectrum) of the kernel k {t, vj). The above-mentioned spectrum, together 
with \i is associated with the eigenfunctions ^^ (t) such that 

1 
^,{t)^X,Jk{t,rj)if^ir])drj. (3.17) 



Numerical realization of the resolvent H{t,ri,X) is a special topic. On the 
whole, there exist a number of efficient algorithms for the solution of Fredholm 
integral equations of the second kind. In any case, both the parameter A and 
the function H(t,ri,X) are very convenient for analytical purposes. There are 
publications devoted to the study of the structure of the resolvent in the neigh- 
borhood of the characteristic numbers. Analytical estimates of the order of the 
growth of the characteristic numbers depending on the properties of the kernel 
are worked out [15] (PP- 61, 62). 

The following example clearly demonstrates the influence of the parameter 
A on the solution. Thus, the equation 

1 



for A 7^ —6 ± 4\/3, has the unique solution [30j (pp. 29, 30) 

1 
fit)^q (t) + ^, ^ ^^^ _ ^^ / [6 (A - 2) {t + ii) - Utrj - 4A] q (r/) dr,. 
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Equation 



1 

^{t) = X [e*'\{r,)drj + q{t), (3.18) 





for A 7^ 1, has the solution 

V{t) = q (t) + Y^e* J e-*q (t) dt- 

if, however, A = 1 and J e^^q (t) dt — 0, then 

^it)=q{t)+ce\ 

where c is an arbitrary constant. For A = Ai, equation (|3.18|) is, in general, 

unsolvable, whereas the corresponding homogeneous equation has an infinite 

number of nontrivial solutions [27 (pp. 74-76). 

As is pointed out by V. I. Smirnov [31] (p. 130), when considering integral 

equations whose kernels are analytical functions of the parameter, one can meet 

with substantial deviations from the regularities present in general theory of 

Fredholm integral equations of the second kind. From this point of view, of 

interest is the paper by Z. I. Khalilov i32j devoted to investigation into the 

equation 

1 

(^ (t) = / [fco (t, 77) + ^lkl (t, ry)] (^ (r;) drj + q (t) , (3.19) 



where /i is a parameter, under certain restrictions on /cq (t, 77), fci {t, 77), and q (t). 

Tamarkin's theorem is given: Equation p.l9p either has no solutions for 
all q (t) and any value of /i, or it has a unique solution for all q (t) and any ^ 
except for, may be, a countable set of values of fj, together with a limit point 
(if it exists) at the infinity. The statement of this theorem for a more general 
formulation of k {t, 77, /j.) is given in [3T] (p. 132). 

In the first case, the author [32] has termed the kernel of equation (|3.19p 
exceptional, whereas in the second case, he has termed it non-exceptional. A 
theorem is proved stating that, in the case of a non-exceptional kernel, equation 
p.l9p can always be reduced to a usual integral equation. Representations of 
the kernel of this equation r (t, 77) and of the free term via the functions and the 
parameter from p.l9p . including the resolvent, are given. 

One more theorem: The spectrum of the kernel in (j3.19p either coincides 
with the whole plane of the complex variable fi or coincides with the spectrum 
of the kernel r {t, rj). Conditions of the solvability of equation (|3.19p in the case 
of multiple eigenvalues are also obtained. 

In the context of this consideration, the material presented by M. L. Krasnov 

[55] is relevant to the case. Thus, it turns out that the homogeneous integral 

equation 

1 

V{t)^ I [p it) p (77) + pa (t) p (77)] ip (77) drj, (3.20) 
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where p (t) and a (t) are continuous functions obeying the relations 

1 
P^ iv) dr] = l; IP iv) o- (??) dr) = 0, 



has the nonzero solution 

^{t) = p (t) + fia (t) 

for any /i [33] (p. 272). In other words, equation p.20p does not represent a 
Fredholm integral equation of the second kind. 

F. Tricomi summarizes: "In the theory of equations that do not obey Fred- 
hom's theorems, three new phenomena often occur: 

- the presence of finite limit points of the spectrum of characteristic numbers 
or even of a continuous spectrum, i.e., of characteristic numbers that fill the 
whole interval of the A-axis and even the whole A-axis; 

- the presence of characteristic numbers of infinite multiplicity, i.e., of char- 
acteristic numbers associated with an infinite number of linearly independent 
functions; 

- the presence of bifurcation points (in the real nonlinear case), i.e., of those 
points of the A-axis by passing through which the number of solutions of the 
equation changes while remaining finite" [57] (pp. 207, 208). 

Note that the above-mentioned nonlinearity is not necessary for the appear- 
ance of various miracles in the behavior of the solutions. As a matter of fact, we 
initially have equation (|3.16p with the function (p (t) that is explicitly present. 
Exactly this type of equations is necessary for the study of peculiarities of the 
solutions, such as branching and bifurcation (they can be interpreted in different 
ways). 

In the general case, to obtain an equation of the second kind, by analogy 
with (|2.16i l one needs to have a nonlinear initial object [34] (pp. 22, 23, 304-305, 
318). By definition [33], the values of A^ from (|3.17p representing characteristic 
numbers of the continuous kernel k (t, rj) are also points of bifurcation of equation 

dsini). 

However, only rather special cases of the dependence of the kernel on the 
parameter have been considered. Thus, Tamarkin's theorem does not envisage 
the presence, alongside with /Lt, of the parameter A. At the same time, if one is 
motivated by the development of an objective basis of analytical investigation, 
it is desirable to make the number of the parameters in the kernel sufficiently 
large. From a formal point of view, there are no obstacles on this way. The 
main question is whether adequate computational and theoretical means exist. 

On the basis of a research into literature, we can arrive at a conclusion that 
during the period whose beginning dates back to the publication [3S] there has 
been no more or less significant progress in the field of the theory of integral 
equations with parameter-dependent kernels. How this fact can be explained? 

It seems that the following answer can be given. On the one hand, there 
are, obviously, considerable difficulties of fundamental character on the way of 
the development of mathematically rigorous theory. On the other hand, the 
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theory of integral equations has become part of functional analysis, and, at the 
modern stage, its problems are considered from a more abstract point of view 
(the theory of an implicit function). In a certain sense, the factor of nonlinearity 
may prove to be of minor importance here [34j (section 1). 

In light of the above, we think that it would be rather interesting to develop 
the theory of Volterra integral equations of the second kind with the parameter- 
dependent kernel. On the one hand, the theory of equations of this type is 
objectively simpler that of Fredholm. On the other hand, such equations are 
associated in a natural way with the procedure of the reduction of Leontief 's 
model to the integral form: see sections 3.1 and 3.2. 

Remaining within the framework of the objective orientation of the present 
consideration, we can draw the following conclusions: 

- the solution to an integral equation of the second kind even in a rather 
particular case of dependence of the kernel on the parameter exhibits a wide 
range of possibilities of behavior; 

- as is illustrated by the simplest examples, the parameter can put the equa- 
tion outside the class of Fredholm equations of the second class that are char- 
acterized by stable dependence of the solution on the data of the problem; 

- such properties of the seemingly linear integral equation of the second kind 
with the parameter-dependent kernel characterize it as a rather efficient tool for 
modelling of dynamic processes; 

- at the same time, the breadth of the class of solutions results in difhculties 
of an analytical study and of numerical realization of mathematical models even 
in the case of dependence on the parameter of a particular type. 
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SUMMARY 



Traditional models of macroeconomic dynamics are fundamentally incorrect. 
The reason lies in a misunderstanding of peculiarities of the analysis of infinites- 
imal quantities. However, even those types of solutions that are envisaged by 
the above-mentioned models are nonrepresentative in the sense of the reflection 
of realities. It became obvious that the techniques of the theory of linear dif- 
ferential equations were insufRcient here. Accordingly, the scientists' attention 
switched to the theory of nonlinear differential equations. 

At the same time, balance and, accordingly, the model with matrix proper- 
ties are objectively inherent in the economic system. For the reduction of this 
model to a differential form, there exist rather elementary means that proved to 
be unclaimed. Macroeconomic rhetoric - the power of the accelerator, a lag on 
the part of demand, etc. - accompanied by the use of a lot of abstract coefficients 
prevailed. 

Why such an entourage? One cannot but get an impression that the issue, in 
essence, is political and is deeply rooted. Economics is seemingly a major science 
because it is more intimately related to the capital than other sciences and it 
works out wise recommendations for the ruling elite. Therefore, economics 
is prescribed to use higher mathematics whose peak is differential equations: 
this has been a mentality since long ago. Nobody has heard anything about 
other equations, whereas operating the categories of balance is inappropriate 
because it looks like school arithmetic. The elite, in its turn, acts according to 
prescriptions of science rather than at will. 

However, there is no organic interrelation between matrix and nonlinear 
differential equations. On the contrary, it can be said that linear theory of 
integral equations originated in matrix analysis. The Fredholm linear integral 
equation of the second kind with a parameter-dependent kernel proves to be 
rather representative with regard to the class of possible solutions. It seems 
that it can be used for the description of any zigzags of the economy. The price 
one has to pay for this is the nontriviality of existing theory. 
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